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ABSTRACT 

There is cumulative evidence showing that, for the most massive galaxies, the fraction 
of disk-like objects compared to those with spheroidal properties increases with red- 
shift. However, this evidence is thus far based on detailed morphological analyses of 
these objects which ultimately rests on the shape of their surface brightness profiles. 
To explore the consistency of this scenario, it is necessary to measure the dynamical 
status of these galaxies. With this aim we have obtained near-infrared integral field 
spectra in the H-band for 10 massive galaxies (AT* ^ 10 n h^Q : M Q ) at z ~ 1.4 with 
SINFONI at the VLT. Our sample is selected by their stellar mass and EW[OII]> 15 
A, to secure their kinematic measurements, but without accounting for any morpho- 
logical or flux criteria a priori. Through this 3D kinematic spectroscopy analysis we 
find that half (i.e. 50±7%) of our galaxies are compatible with being rotationally sup- 
ported disks, in agreement with previous morphological expectations. This is a factor 
of approximately two higher than what is observed in the present Universe for objects 
of the same stellar mass. Strikingly, the majority of our sample of massive galaxies 
show clear and fairly large rotational velocity maps, implying that massive galaxies ac- 
quire rapidly rotational support and hence gravitational equilibrium. Our sample also 
show evidence for ongoing interactions and mergers. This demonstrates that massive 
galaxies at this epoch are not a homogenous sample, but show significant diversity and 
must have continued evolution beyond the fading of stellar populations, to become the 
massive galaxies we find in today's Universe. 

Key words: galaxies: evolution - galaxies: massive galaxies - galaxies: high-redshift 
- galaxies: 3D spectroscopy 



1 INTRODUCTION 

Massive (M ste iiar W 11 h^g Mq) galaxies represent a chal- 
lenge to the dominant ACDM paradigm as their proper- 
ties (such as number densities, star formation histories, size 
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growth, etc.) should be well constrained by galaxy evolu- 
tion models, but seemingly this is not always the case (e.g., 
Benson et al. 2003, Baugh et al. 2006; Conselice et al. 2007, 
Trujillo 2012). Several recent studies find that these galax- 
ies are better represented by a disk-like population at z=2-3 
(Cameron et al. 2011, Van der Wei et al. 2011, Bruce et al. 
2012, Buitrago et al. 2013) as opposed to their counterparts 
in the local Universe, which are preferentially large early- 
type galaxies (Baldry et al. 2004). A dynamical test of this 
scenario is still missing, and would complement the great 
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deal of morphological information gathered to date (Oesch 
et al. 2010, Weinzirl et al. 2011, Wang et al. 2012, Van der 
Wei et al. 2013 in prep.). 

Observationally, most of the information regarding high 
redshift galaxies comes from large and deep NIR surveys, 
which aim to probe galaxy evolution during the last 12 
Gyr. The data gathered so far favours a picture in which 
late- type and clumpy /interacting objects are more common 
at high redshift, due to the higher gas fractions for these 
systems (Erb et al. 2006, Tacconi et al. 2010) and higher 
merging rates (e.g. Conselice et al. 2008, Lopez-Sanjuan et 
al. 2011, Bluck et al. 2012). Ultimately, we would ideally 
like to rely on spectroscopic information to fully test and 
characterize all the processes involved in galaxy assembly. 
However, even for massive (and thus very often luminous) 
galaxies, it is very expensive in observational time to obtain 
high signal-to-noise ratio spectra. From the tens of massive 
galaxies studied so far at z > 1.5 with traditional long-slit 
techniques (Kriek et al. 2006, Cimatti et al. 2008, Newman 
et al. 2010, Onodera et al. 2010, van de Sande et al. 2011, 
Toft et al. 2012) there is some agreement on the high velocity 
dispersion values of these objects, confirming their inherent 
massive nature even at such early cosmic times (Cenarro & 
Trujillo 2009, Cappellari et al. 2009, Van de Sande et al. 
2012). Despite this, the true nature of these systems and 
how they evolve is not yet understood. 

Integral field spectroscopy (aka 3D spectroscopy) is 
presently a well-established technique which can enhance 
greatly our understanding of these massive galaxies, both 
at low-redshift (Cappellari et al. 2011, Sanchez et al. 2012) 
and at high-redshift (Law et al. 2009, Forster-Schreiber et 
al. 2009, Contini et al. 2012). On the one hand, it measures 
the rotational and the velocity dispersion support (currently 
from gaseous kinematics) for a galaxy, and thus provides 
us with physical information, such as both baryonic and 
dark matter components, as well as the presence of rotation. 
Spatially-distributed spectral information furthermore helps 
us address the question of how morphology and galaxy as- 
sembly are linked. State-of-the-art Integral Field Unit (IFU) 
studies reveal a kinematic mixture within high-z systems 
(e.g., Genzel et al. 2008, Nesvadba et al. 2008, Shapiro et al. 
2008, Cresci et al. 2009, Law et al. 2007, 2009, Wright et al. 
2007, 2009, Epinat et al. 2009, Bouche et al. 2010, Lemoine- 
Busserolle & Lamareille 2010, Gnerucci et al. 2011, Epinat 
et al. 2012), with a high percentage of objects displaying 
large ordered rotational motions, and also fairly large veloc- 
ity dispersions. 

In the low redshift Universe, the SAURON survey (de 
Zeeuw et al. 2002) and the more recent ATLAS 3D survey 
(Cappellari et al. 2010) have opened a new perspective on 
the kinematics of local massive galaxies. They classify early- 
type galaxies as slow and fast rotators depending on the de- 
gree of angular momentum they exhibit. Fast rotators have 
been found to host disks made up of gas and stars which con- 
tain a range of the galactic mass fraction (Falcon-Barroso et 
al. 2006, Krajnovic et al. 2008, Davis et al. 2011, Krajnovic 
et al. 2012). In the other hand, at high-z, not only mas- 
sive galaxies seem to be better represented by a disk-like 
population (e.g., Buitrago et al. 2013) but, for those with 
early-type morphology, a flatter nature (possibly indicating 
a disk-component) is prevalent (Van der Wei et al. 2011, 
Chevance et al. 2012, Chang et al. 2013a, Chang et al. 2013b 



in prep.). Our study aims to provide data on the kinematics 
of massive galaxies at intermediate redshifts (z ~ 1.4) and 
thus covering the gap between current and distant cosmic 
epochs. 

We discuss in this paper IFU studies of a sample of 
galaxies with log M„ > 11 at high (z ~ 1.4) redshift. These 
are good targets for IFU spectroscopy due to their relative 
compactness (which make them easy to observe even with a 
small field-of-view) and the current lack of spectra for this 
galaxy population. Ideally, absorption lines measurements 
would be the best indicators to examine the stellar popula- 
tions and the motions within this population (e. g. Bedre- 
gal et al. 2009). Nevertheless, Ha emission line analyses are 
usually preferred due to its relative ease of study. It is very 
important to emphasize the fact that one must be cautious 
when interpreting the derived Ha kinematics, since this ion- 
ized gas is collisional and dissipative and it may not be cou- 
pled with the stellar component in the galaxy. However, it 
is not unreasonable to expect good correlation between gas 
emission and broadband imaging (e.g. Forster-Schreiber et 
al. 2011), especially when studying relaxed systems. 

We present in this paper 8m-telescope VLT observa- 
tions of these massive galaxies using the SINFONI IFU spec- 
troscopy instrument. Our sample comprises 10 objects at 
z ~ 1.4, whose redshift choice is a trade-off between high 
redshift and having a sample selected by stellar mass. Nev- 
ertheless, this redshift is of particular importance in the de- 
velopment of massive galaxies in particular, and galaxies in 
general, as it nearly coincides with the star formation and 
quasar activity peak, apart from being the epoch at which 
massive galaxies begin to switch their morphologies from 
late to early types (Buitrago et al. 2013). We discuss the Ha 
derived kinematic properties for this sample and interpret 
these findings in the context of massive galaxy formation. 

This paper is structured as follows: Section [2] describes 
the data reduction and subsequent analysis, Section [3] dis- 
cusses the different tests 3D spectroscopy offer for character- 
izing the rotation dominated nature of these objects, and in 
Section |4] we present the conclusions of our study. There is 
a final Appendix in which we show each galaxy, explaining 
in a detailed way its particularities, and how every mas- 
sive galaxy is related with the total sample. We name the 
galaxies in our sample with the prefix POWIR followed by a 
number, instead of the numeric code in the parent survey, as 
it simplifies tracking the individual galaxies throughout the 
paper (please note there is not a POWIR9 galaxy). Besides, 
these numbers are plotted along with the galaxy symbols 
aiming at the same purpose. We assume a concordant cos- 
mology: Ho=70 km s~ Mpc -1 , Q\ = 0.7, and S7 m = 0.3, and 
use a Chabrier (2003) IMF and AB magnitude units unless 
otherwise stated. 



2 DATA AND ANALYSIS 
2.1 Observations 

The parent sample from which our target galaxies are se- 
lected comes from the Palomar Observatory Wide InfraRed 
survey (POWIR; Bundy et al. 2006, Conselice et al. 2007, 
2008). This survey covers a 1.53 deg 2 area in the K s and J 
bands down to Kvega = 21 and Jv ega = 23.5. This imaging 
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Table 1. Observational data for our massive galaxy sample 



Name 


POWIR ID 


RA(J2000) 


DEC(J2000) 


z 


log M* 


K s mag 


Observ. night 


Integration time 


(S/N),„ ops 


Seeing 












log (h^M 3 ) ' 


/ega magnitudes 




sec 






(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


POWIRl 


32007614 


23.516207 


0.043184 


1.374 


11.2 


18.2 


21 — Jun 


5400 


3 


0.58 


POWIR2 


32073051 


23.519257 


0.111119 


1.396 


11.0 


18.6 


04 - Sep 


5400 


3 


0.52 


POWIR3 


32015443 


23.500269 


0.146335 


1.384 


11.2 


18.5 


08 - Aug 


2700 


3 


0.65 


POWIR4 


32015501 


23.501966 


0.102931 


1.394 


11.4 


18.0 


31 - Jul 


2700 


2 


0.58 


POWIR5 


32021317 


23.499363 


0.161825 


1.382 


11.3 


18.3 


20 - Jul 


5400 


2 


0.55 


POWIR6 


32021394 


23.496317 


0.159318 


1.375 


11.5 


17.9 


01 - Jul 


5400 


3 


0.48 


POWIR7 


32029850 


23.517859 


0.282879 


1.396 


11.3 


18.3 


19 - Jul 


5400 


2 


0.63 


POWIR8 


32037003 


23.503817 


0.334799 


1.400 


11.0 


18.4 


01 - Jul 


5400 


3 


0.42 


POWIR9 






















POWIR10 


32100059 


23.516451 


0.024167 


1.389 


11.1 


18.7 


20 - Jul 


5400 


3 


0.55 


POWIRl 1 


32100778 


23.510761 


0.235435 


1.393 


11.1 


18.3 


08 - Aug 


2700 


2 


0.65 



Notes. (1) Name of the galaxy (2) Name of the galaxy in the parent POWIR/DEEP2 survey (3) Right ascension (4) Declination (5) 
Spectroscopic redshift from our SINFONI observations (6) Stellar mass from the parent POWIR/DEEP2 survey (7) ii" s -band magnitude 
from the parent POWIR/DEEP2 survey (8) Date of the observations, in 2007 (9) Integration time (10) Signal-to-noise ratio threshold 
above which we show the spaxels in the kinematic maps (11) Seeing as derived from the SINFONI telluric standards observed. 



consists of 75 Palomar WIRC camera pointings, with a pixel 
scale of 0.25" /pix. In the K s band 30-s exposures were taken, 
for total l-2h integration time per pointing, and with typical 
seeing of 0.7— 1". In addition, optical coverage was supplied 
with the 3.6m Canada-France-Hawaii Telescope (CFHT) us- 
ing the CFH12K camera in the B, R and I bands. The R- 
band depth is Rab = 25.1, with similar results for the other 
two bands (see Coil et al. 2004a for more details). Both 
Palomar and CFHT images were analysed using 2" diame- 
ter apertures. 

Stellar masses were derived with the photometric tech- 
niques discussed in Bundy et al. (2006) using a Chabrier 
(2003) IMF. Our stellar mass computational method con- 
sists of fitting a grid of model SEDs constructed from 
Bruzual & Chariot (2003) stellar populations synthesis mod- 
els. The star formation histories are parametrized by utiliz- 
ing so-called tau-models (SFR oc exp(— t/r); with r ran- 
domly selected from a range between 0.01 and 10 Gyr, and 
the age of the onset of star formation ranging from to 10 
Gyr), with a range of metallicities (from 0.0001 to 0.05) and 
dust contents (parametrized by the effective V-band optical 
depth tv, where the used values were ry = 0.0, 0.5, 1, 2). To 
analyze the impact of TP-AGB emission, the same exercise 
was also performed with Chariot & Bruzual (2007) models, 
inferring slightly smaller masses (on the order of ~ 10%). 
It turns out that parameters such as metallicity, e-folding 
time or age are not as well constrained as stellar mass due 
to the various degeneracies. The final error in stellar mass is 
measured as 0.2-0.3 dex, i.e. roughly a factor of two (Bundy 
et al. 2006, Conselice et al. 2007, Griitzbauch et al. 2011). 

Spectroscopic redshifts were measured within the 
DEEP2 Redshift Survey (Coil et al. 2004b) using the 
DEIMOS spectrograph (Faber et al. 2003) at the Keck 
II telescope. Spectra were obtained with a resolution of 
R ~ 5000 within the wavelength range 6500 - 9100 A. Red- 
shifts were measured comparing templates to the data, and 
we utilised only those in which two or more lines were iden- 
tified. We show in Figure [T] how our massive galaxy sample 
compares to the rest of the DEEP2 spectroscopic sample. 
It is noteworthy that our final choice of galaxies to observe 
was not based on a selection using colours, morphologies 
or sub-mm flux as many other IFU studies. Additionally, 
we have not only explored compact massive galaxies. In- 
stead our sample is solely selected by stellar mass without 
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Figure 1. Mass-redshift relation for the DEEP2 spectroscopic 
sample. Our massive galaxies are highlighted by the red squares. 
Our study is an attempt to characterise the high redshift high 
mass population of this spectroscopic survey. 



any other criteria a priori except with an [O II] equivalent 
width > 15 A. This last condition might bias our sample 
towards investigating star forming systems. However, this is 
a necessary requirement to robustly asses our galaxy sam- 
ple's kinematic features. We must stress that most massive 
galaxies found at similar redshifts than our sample (z ~ 1.4) 
have significant star formation rates, as highlighted by var- 
ious studies (e.g., Conselice et al. 2007, Perez-Gonzalez et 
al. 2008, Bauer et al. 2011). Adding this information along 
with the fact that various morphologies are found in this 
work suggests that we are close to having a representative 
sample selected by stellar mass. 

Our group was granted 20 hours of observing time in 
service mode with the 3D-spectrograph SINFONI (Eisen- 
hauer et al. 2003, Bonnet al. 2004) at ESO-VLT located 
in UT4-Yepun. Our observations were conducted during 9 
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nights from June to September 2007 - ESO run ID 079.B- 
0430(A) -. SINFONI was used in seeing limited mode and 
thus with a spatial sampling of 0.125" x 0.25". Due to the 
redshifts of our sample of galaxies we choose to observe them 
in the H-band in order to map the Ha emission. We also 
took special care in our final galaxy selection to ensure that 
their Ha emission lines were not situated close to any OH 
sky emission lines - based on the atlas from Rousselot et al. 
(2000) - which would potentially hamper our results. The 
spectral resolution (R ~ 3000) allows us to disentangle sky 
emission lines close to our target. 

Our observational strategy was the so-called 'butter- 
fly pattern' or 'on-source dithering', by which the galaxy is 
set in two opposite corners of the detector to remove sky 
background using contiguous frames in time. Several galax- 
ies in our sample (POWIR4, POWIR5 and POWIR7) could 
only be observed half of their nominal integration time (lh 
30min). Even in these cases, exquisite SINFONI sensitivity 
permitted us to detect the emission from all our objects. Im- 
ages were dithered by 0.3" in order to minimize instrumen- 
tal artefacts when the individual observations were aligned 
and combined together. PSF and telluric stars were also ob- 
served along with each galaxy for calibration purposes. Mea- 
sured PSFs are listed in Table [TJ for a mean seeing of 0.56" 
throughout our observations. 



2.2 Data reduction & observed kinematic maps 

We have used the ESO-SINFONI pipeline version 2.5.0 
(Modigliani et al. 2007, Mirny et al. 2010) to reduce our 
data. In brief, this pipeline subtracts sky emission lines 
(using algorithms by Davies et al. 2007), corrects the im- 
age using darks and flat-fields, spectrally calibrates each 
individual observation and reconstructs all the information 
into a final datacube. The recipe used for this purpose was 
sinf o_rec_j itter, which was fed exclusively with the mas- 
ter files provided by ESO. All of these processes were per- 
formed separately for each individual exposure. Afterwards 
the two datacubes were combined into a single one by us- 
ing the recipe sinf o_utl_cube_combine. We always used the 
pipeline parameter product-density — 3 (which retrieves the 
most detailed possible outputs), objnod-scalessky — true (to 
perform a subtraction of the median value at each wave- 
length and thus remove the sky more efficiently) and sky- 
cor.rot_cor = true (to remove the contribution of any rota- 
tional OH transitions). 

The final datacube was spatially smoothed using a sub- 
seeing Gaussian core (FWHM=2 pixels) to increase the 
Signal-to-Noise Ratio (SNR) without affecting our data in- 
terpretation. We analysed this datacube with IDL routines 
we constructed. Basically, we located the Ha line in each 
spaxel according to the known spectroscopic redshift of the 
target galaxy, and then fit a Gaussian profile, taking into 
account the sky spectrum weighting its contribution with 
the help of the routines mpfit and mpfitfun (Markwardt 
et al. 2008). Radial velocity maps were computed using the 
relativistic velocity addition law: 

(^spaxel ^cen) 

spaxel — z j C 

J- r Zcen 



where z spaxe i and z cen are the redshifts for a given spaxel 
and for the kinematic centre of the galaxy, respectively. 

From the Ha line width, we computed velocity disper- 
sion maps, subtracting the instrumental broadening, mea- 
sured from sky lines. In addition, we obtained Ha and [Nil] 
A6583A line flux maps. The SNR per spaxel was calculated 
in the following manner: the signal was the intensity of the 
Ha line, and the noise was the standard deviation of the 
residual spectrum, with both signal and noise weighted by 
the sky contribution around the Ha wavelength. A contin- 
uum map was also constructed with the spectral information 
in the range 1.5 — 1.7/xm, i.e., in all the H-band except its 
borders, where the information was noisier. We fit a linear 
function to the galaxy spectrum in this wavelength range 
in order to account for the existence of continuum emission 
and its possible variation within this wavelength range. Our 
continuum maps show the integral of the fitted mathemati- 
cal function. Spaxels with values below zero are coloured in 
white. 

All our objects have continuum emission which we have 
compared to the ionized gas emission, which in principle tells 
us only about the areas of star formation in each galaxy. As 
such, comparing the peak of the emission and the continuum 
gives us insights to how well Ha , and hence star formation, 
traces the underlying older stellar populations. One caveat 
to this was a known problem with the SINFONI detector, 
where there are stripes of flux in the data after coadding high 
numbers of spectral pixels. One can see these stripes for in- 
stance in the continuum maps of POWIR1 (see Appendix 
A; southern part of the galaxy) and POWIR3 (see Appendix 
A; white stripe on the top of the galaxy). This effect pre- 
vented us from making a total continuum flux measurement, 
but did not affect qualitatively the fact that we could locate 
where the maximum of the continuum was in the detector. 

2.3 Data modelling 

We recovered the kinematic parameters for each galaxy in 
our sample by fitting a model to the velocity field obtained 
from our SINFONI datacubes. To perform this task we as- 
sumed that the galaxies from our sample are described kine- 
matically as rotating disk systems with a symmetric rotation 
curve. For this method we utilised the formalism and pro- 
grams developed in Epinat et al. (2009). The full theoretical 
description of this method is in Epinat et al. (2010), where 
the authors also conducted a comparison with local galax- 
ies to asses the reliability of their method. Essentially this 
consists of a \ 2 minimization between the observed data 
and a given high resolution model convolved to our pixel 
scale and seeing conditions. We chose the flat rotation curve 
parametrization used by Wright et al. (2007, 2009), as sug- 
gested in Epinat et al. (2010) from the study of local galaxy 
velocity fields projected at high redshift. The analytical ex- 
pression for this is given by: 

V(r)=V t - 
r t 

when r ^ r t and 

V(r) = V t 

otherwise. In the above equations Vt is the value for the 
plateau in the rotation curve and r t is the radius at which 
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the plateau is reached. The model contains seven parameters 
: the center (x c and y c ), the systemic redshift (or velocity), 
the inclination of the disk, the position angle of the major 
axis and the two rotation curve parameters. Note that the 
fit to these simple formulas were done by considering the 
associated error map for the velocity field. 

As shown and discussed in Epinat et al. (2010), due 
to the reduced spatial information of our data, and due to 
some degeneracy in the models, the center and the inclina- 
tion are the least constrained parameters. We thus fixed the 
center to the spaxel with the maximum flux in the contin- 
uum maps (in principle the continuum may trace better the 
galactic center) as well as the inclination, reducing to four 
the number of free parameters in our model. In rotating disk 
models there is also a degeneracy between rotation velocity 
and inclination (its sine) that can only be solved using very 
high resolution data. As a result the inclination is the major 
source of uncertainty for determining the actual rotation ve- 
locity. Given the photometric quality of our POWIR parent 
sample imaging, it is difficult to constrain this parameter 
with a high certainty. We used GALFIT (Peng et al. 2010) 
to fit single Sersic (1968) surface brightness profiles to our 
sample. We have only utilised the inclination retrieved by 
the GALFIT analysis, and discarded the rest of the struc- 
tural parameters, owing to the high uncertainties given the 
large Point Spread Function (PSF) of our images. We used 
bright, non-saturated stars within our imaging as a model 
PSF. The output inclination was then utilised as input for 
our velocity modelling. 

Once the best model for the rotational velocity was ob- 
tained, we also computed a model velocity dispersion map. 
To calculate this we took into account the width of the 
Ha line due to the unresolved velocity gradient. The intrinsic 
velocity dispersion is obtained after subtracting in quadra- 
ture the velocity dispersion model map from the observed 
one. Finally, to facilitate a comparison with other samples 
and to discuss each galaxy as a whole, we also computed 
the velocity dispersion value doing a weighted average of the 
value of every spaxel by an amount inversely proportional 
to their squared error. Note also that this velocity disper- 
sion is computed from the beam smearing corrected maps. 
To avoid confusion, henceforth we will refer to this quantity 
as the 1/error 2 velocity dispersion. Results from the models 
are listed in Tableful and their maps are the first two mon- 
tages in the Appendix A. All of the galaxy kinematic maps 
are also in the appendix. 



3 RESULTS 

3.1 Kinematic classification 

Previous 3D spectrocopy studies of high redshift galax- 
ies have demonstrated that these systems are more 
clumpy /irregular and have higher velocity dispersions than 
local galaxies (Forster-Schreiber et al. 2011 and references 
therein). There are several attempts in the literature to es- 
tablish a kinematic classification of high-z galaxies relying 
on Ha kinematics (Flores et al. 2006, Law et al. 2009, Cresci 
et al. 2009, Forster-Schreiber et al. 2009, Epinat et al. 2009, 
Gnerucci et al. 2011, Epinat et al. 2012). All of these studies 
roughly agree that there are three basic kinematic classes, 



which may be linked with the morphological nature of each 
galaxy. First, rotating disks have been observed, showing 
well-defined and regular rotational velocity gradients that 
are larger than their velocity dispersion. Usually these sys- 
tems are large in size. Ongoing mergers are also clearly dis- 
tinguished, not only by disentangling two or more compo- 
nents but through a chaotic velocity pattern, and local incre- 
ments in the velocity dispersion. Finally, objects which do 
not fit in any of the previous categories are tagged as per- 
turbed rotators, which are probably more similar to early 
type systems due to their high velocity dispersion in com- 
parison with their maximum rotational velocity. 

Before characterising our sample according to these cri- 
teria, we enumerate several caveats that might affect our 
interpretation of the data. First, behind this classification 
there is the disk-like assumption in the modelling. This will 
not be an accurate model when dealing with mergers or 
pure spheroidal galaxies. Arguably, this has an impact in 
our sample as massive galaxies in the local Universe are 
predominantly early types (e. g. Baldry et al. 2004). How- 
ever, at the redshift of our observations (z ~ 1.4), we would 
expect to have a morphological mixture (van Dokkum et al. 
2011, Buitrago et al. 2013). Finally, we must remember we 
are looking at the gas emission and not at the total stellar 
component. 

Performing a preliminary investigation, we notice that 
the Ha emission extends over all of the K s and continuum 
images for most cases, which is difficult to reconcile with 
the possibility of being spheroid-like objects, specially when 
Ha emission is usually linked with star formation. As stated 
in the introduction, the existence of gas disks within early 
type galaxies has been reported in the past (Falcon-Barroso 
et al. 2006, Krajnovic et al. 2008, Oosterloo et al. 2010, 
Krajnovic et al. 2012) but their sizes (hundreds of parsecs) 
are much smaller than our current gas disks which span the 
galactic size. In fact, these disks are a fundamental aspect 
for understanding the evolution of these galaxies. 

We focus now on the individual properties of the sam- 
ple, which has been discussed at length in the Appendix A. 
In summary, POWIR1, 2, 6 and 8 show clear rotational gra- 
dients. However, the centers of the latter two display fairly 
large velocity dispersions. In both Epinat et al. (2009) and 
Forster-Schreiber et al. (2009) these were indications of disks 
galaxies developing a bulge component. It is also notewor- 
thy to mention that, in case of a very steep gradient or ring 
in Ha , a central peak due to beam smearing is expected. 
However, these gas distributions are mostly found in the 
local Universe in galaxies with large bulges (Peletier et al. 
2007), and thus this could be in agreement with a progres- 
sive evolution towards early-type systems. Nevertheless, in 
light of the kinematic probes, we classify the four galaxies 
as rotating disks. 

On the other hand, POWIR4 and 10 are classified as in- 
teracting objects. In fact, in both cases the Ha emission orig- 
inates (although there is a very weak detection in POWIR10 
for the main object) from spaxels that do not belong to 
the target galaxy. Hence we discern two separate galax- 
ies interacting within our SINFONI data. The photometric 
data has been derived for the main objects, which we iden- 
tify as the massive objects in the K s band imaging, while 
the Ha detection comes from the secondary galaxies. These 
non massive galaxies are therefore excluded from our dis- 
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) galaxies in our sample, also called the anisotropy plot. Apart from 



these, we supplemented the figure with the local early-type galaxies in Emsellem et al. (2011). These low redshift objects are early-type 
galaxies studied as part of the ATLAS 313 survey (Cappellari et al. 2011). Ellipticities for our sample were measured in the K-band 
imaging of POWIR/DEEP2 survey using GALFIT and thus taking into account the PSF of our imaging. The continuous line defines 
the ideal oblate rotator with isotropic stellar velocity distribution for integral field studies (Binney et al. 2005, Cappellari et al. 2007). 
The uncertainties are large, but it is clear that massive galaxies at z ^ 1.4 depart from the low-z velocity dispersion dominated objects. 



cussion and plots, but we derived kinematics for them to 
understand which physical processes are taking place in 
the merger. Little can be said about the two massive and 
main objects: POWIR4 is completely devoid of Ha emission, 
while POWIR10 looks like a point source with a strong 
[Nil] line in its center, that suggests it hosts an AGN. 
POWIR3 is an object which we include in this category as 
well, as its K s and continuum images do not match with the 
Ha emission, whose map is quite irregular. 

The rest of the objects are more difficult to catalog visu- 
ally. We must bear in mind that 2 out of the 3 other objects 
were observed half of the nominal integration time. Either 
POWIR5, 7 and 11 have relaxed morphologies in the K s and 
continuum bands while Ha shows, as expected, a more com- 
plicated pattern. POWIR11 is different, despite the obser- 
vational issues. It has an easily distinguishable and large ro- 
tational field, which fits better the disk object classification. 
The other two galaxies are catalogued as perturbed rotators. 
One of our major conclusions is that these massive systems 
have a wide range of properties when examined kinemati- 



cally, and that many of them display significant rotational 
support. 

3.2 Rotation vs velocity dispersion dominance 

In low redshift studies, the V m ax/<J vs. e diagram (also called 
the anisotropy plot; Figure [2} is a classical tool to measure 
the kinematics of early type galaxies (Illingworth 1977, Ben- 
der et al. 1994, Cappellari et al. 2007, Emsellem et al. 2011). 
We created this plot with our sample's data as an exercise, as 
massive galaxies at low-z are predominantly early type sys- 
tems and therefore this is a good test to shed light into the 
nature of our sample. However, we remind the reader that 
we are dealing with ionized gas emission instead of stellar 
kinematics. 

The plotted parameters used in this relation usually are 
measured at one effective radius distance from the galaxy 
center. To palliate our uncertainty on this, we computed 
effective radii in our sample using the relation published 
in Buitrago et al. (2008) for massive disk-like galaxies (to 
be consistent with our modelling), extrapolated to each 
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Figure 3. Inclination vs rotational velocity, 1/error 2 velocity dispersion and Vmax/a for each galaxy (top row) and for each galaxy 
within one effective radius (bottom row). Galaxies are represented by their respective numbers (for more information see Table[2), having 
the following color coding: blue for disk-like galaxies, red for perturbed rotators and orange for interacting galaxies. Pearson correlation 
coefficients for each galaxy appear in the bottom right corners. In light of these plots, the fact that we do not find any correlation 
indicates there is likely no strong bias on obtaining the kinematic data for our sample based on inclination effects. 



galaxy's redshift. Then we computed Vmax/& in the clos- 
est aperture to the calculated effective radius, as to measure 
velocity dispersions we need an integer number of spaxels 
around our kinematic center. All the information used is 
tabulated in Table [3] We add low redshift galaxies from the 
ATLAS 313 Survey (Emsellem et al. 2011). Note however that 
their kinematics are obtained for the stellar component and 
that not all their masses fulfill our definition of a massive 
galaxy (M» > 10 11 Mq). Nevertheless the comparison makes 
sense as their sample is composed of some of the most mas- 
sive galaxies in the nearby Universe. Although uncertainties 
are large (also for the ellipticity, due to the coarse resolution 
of our images), we find that all of the massive galaxies at 
high-z from our sample lay above the line defined by 




which is the minimal rotational approximation to the 
isotropy line optimized for integral-field kinematic obser- 
vations (Binney et al. 2005, Cappellari et al. 2007). This 
reveals the high level of rotational support for these mas- 
sive galaxies at z~1.5, especially when comparing with slow 



rotators that are the most massive galaxies nearby. It is in- 
teresting as well that both POWIR2 and POWIR8 are close 
to the isotropy line. Both galaxies show a clear disk with 
a velocity dispersion enhancement in the center, which is 
potentially related to a forming bulge component. Possibly 
these galaxies are beginning to transform into an early-type 
classification due to secular and passive evolution. However 
due to their total Vmax/cr (those not derived within one ef- 
fective radius) and rotational velocity fields we acknowledge 
them as more similar to disk dominated galaxies. 

In order to assess whether these findings are biased, we 
constrain the impact of the galaxy inclination in our sam- 
ple. Figure [3] shows how our kinematic velocities (namely 
maximum rotational velocity, 1/error 2 velocity dispersion, 
their respective ratio, and these values within one effective 
radius) are correlated according to their inclinations. Pear- 
son correlation coefficients are supplied in each subplot. The 
small values of the Pearson coefficients confirm they do not 
depend on the galaxy inclination, and thus excluding this 
factor as a dominant source of error. 

Comparing the rotational velocity versus the velocity 
dispersion within our systems is a good way to gauge the 
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Figure 4. Maximum rotational velocity inferred from our modelling versus the 1/error 2 velocity dispersion after correcting it for beam 
smearing. Numbers depict each one of the massive galaxies from our sample, whereas the violet triangles come from the MASSIV sample 
(Epinat et al. 2012) and the green squares from the SINS sample (Forster-Schreiber et al. 2009). Note that, for these latter objects, error 
bars are asymmetric. We also attach the histogram of the Vmax/cr of our massive galaxies with and without the addition of the massive 
galaxies in other samples (dashed or solid histogram respectively). For all these massive galaxies we find that V m ax/& > 1, as they lay 
above the 1:1 solid line, with most of them showing ratios 3—5 which corroborates their gravitational support. The fact that SINS and 
MASSIV objects lay in the upper part of this plot is further evidence that these systems are more disk-like. 



dynamical status of massive galaxies at high redshift, where 
the information is not so detailed as in the local Universe. 
This is a completely separate method for obtaining informa- 
tion on the nature of these galaxies beyond imaging. This 
comparison is shown in FigureUwhere we use the maximum 
rotational velocity from our models, and the 1/error 2 veloc- 
ity dispersion for our sample. We supplemented this infor- 
mation with SINFONI data from published massive galax- 
ies with modelling information available, from the MAS- 
SIV survey (galaxies VVDS140258511 and VVDS220584167 
from Contini et al. 2012, Epinat et al. 2012) and the SINS 
samples (galaxies Q2343-BX610, D3a-6004,D3a-6397, D3a- 
15504 from Forster-Schreiber et al. 2009). 

As can be seen, all the galaxies in these samples exhibit 
Vmax/& > 1, in most cases with values larger than 3. We 
construct as well the histogram of the data shown in Figure 
[4] The dashed part corresponds to the galaxies that are not 
part of our sample. Although the number statistics are poor, 
all the massive galaxies plotted show rotational velocities ex- 
ceeding their computed central velocity dispersions, in most 



cases by a large factor. Interestingly, the objects from the 
SINS and MASSIV surveys have V ma x/cr ratios which are 
on average larger than our values. One possible explanation 
is that, as these objects are selected by their star-formation, 
they are potentially even more rotationally dominated than 
our sample. This is another piece of evidence that several 
massive galaxies in our sample have settled down by z ~ 1.4, 
and are developing a possible bulge component, as suggested 
by the anisotry plot. 

3.3 Dynamical masses 

Our integral field spectroscopy results may also be used to 
explore the dark matter content in our sample. To achieve 
this aim we computed dynamical masses combining the in- 
formation coming from the rotational velocity and the ve- 
locity dispersion maps using the formula (from Epinat et al. 
2009) 

M dyn = Me + M a = + 
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Figure 5. Comparison between the inferred dynamical mass and the photometric stellar mass of our galaxies. Numbers depict each one 
of the massive galaxies from our sample, whereas the violet triangles come from the MASSIV sample (Epinat et al. 2012) and the green 
squares from the SINS sample (Forster-Schreiber et al. 2009). The dynamical mass is obtained using the enclosed mass estimation due to 
the rotational velocity and adding the contribution of the velocity dispersion using the asymmetric drift correction (Meurer et al. 1996). 
However, for the four objects corresponding to the SINS survey, dynamical masses are inferred assuming the isotropic virial estimation 
(see section 9.6 in Forster-Schreiber et al. 2009). The solid line is the 1:1 relationship, while the dashed line is the local relationship for 
local SDSS galaxies in Van de Sande et al. (2011). Assuming good (albeit with 0.2-0.3 dex errors) stellar mass calculations, the difference 
in dynamical mass may relate with the fact that what we measure is the inner ionized gas kinematics in our galaxies and not the overall 
baryonic matter contribution. 



where Me describes the mass enclosed up to a radius Riast 
and M a is called the asymmetric drift correction (Meurer et 
al. 1996) by which we take into account the 1/error 2 velocity 
dispersion (a), i.e., the random motions within the galaxy. 
All of these terms are described in Table [4] except h, which 
is the gas surface density disk scale length described by a 
Gaussian function, which can be expressed analytically as 
h = (21n2)~ 1//2 r c . We assume that both Mg and r e are for 
disk-like systems, as explained in Section[3]2] The outcomes 
of this calculation are plotted in Figure [5] Two lines are 
drawn on this figure: the solid one is the 1:1 reference, while 
the dashed line is the local relationship (Mdyn/M, = 1.68, 
the average ratio for SDSS galaxies) from Van de Sande et 
al. (2011). 

In principle, one would expect all galaxies to populate 
the region above the solid 1:1 line, as their dynamical mass 
would have to account for the baryonic mass plus the dark 
matter component. This is not the case for all the objects 



in our sample. There are a number of reasons which may 
explain this disagreement. First, we must not forget that 
this dynamic mass originates from the ionised gas dynam- 
ics which may depart from the values obtained from the 
stellar measurements. Secondly, our calculations account for 
the mass within Riast, i.e., the maximum Ha radius, which 
is smaller than the apertures where the stellar masses have 
been measured, and also smaller than the typical radius used 
in other works such as Forster-Schreiber et al. (2009) to ob- 
tain this parameter. Our dynamical masses would be larger 
if we correct them for these effects. Furthermore, Martmez- 
Manso et al. (2011) claim a possible overestimation on the 
stellar masses we are utilising from Bundy et al. (2006) and 
Conselice et al. (2007). If confirmed, our dynamical masses 
would be in better agreement with the new stellar masses. 
There could also be present inclination or beam smearing 
correction problems. To try to understand better the origin 
of these dynamical masses, we performed a plot of the rela- 
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Figure 6. Relative importance of the two terms which contribute to the dynamical mass. The number of each galaxy is plotted, to 
better understand the properties of every galaxy. The solid line shows the 1:1 relation for convenience. Note that the MASSIV massive 
galaxies (violet triangles, Epinat et al. 2012, changing its masses accordingly to our Chabrier IMF) are overplotted, but not the SINS 
objects (Forster-Schreiber et al. 2009), as their dynamical mass were calculated assuming the isotropic virial estimation (see section 9.6 
in Forster-Schreiber et al. 2009). The various kinds of galaxies within the present work are spread over this plot, with no obvious split 
in their observational properties. 



tive contributions of each term in the formula used in Figure 
[6] Apart from the low SNR observations of POWIR7, the 
contribution for the velocity dispersion term is quite impor- 
tant for all our sample, making its addition mandatory to 
retrieve correct dynamical mass estimates of high-z galaxies. 

3.4 Tully-Fisher relation 

The Tully-Fisher Relation (TFR) links the maximum rota- 
tional velocity of spiral galaxies with their luminosity or stel- 
lar mass (Tully & Fisher 1977; Fernandez-Lorenzo et al. 2009 
for a comprehensive update). It has also been extended to 
SO galaxies (e.g. Bedregal et al. 2006) and early-type galax- 
ies in general (e.g. Davies et al. 2011). Modern investiga- 
tions focus on finding and understanding any evolution in 
its slope, intercept, or both. It is a powerful scaling rela- 
tion which accounts for how the stellar mass and the dark 
matter content are related (e.g. Conselice et al. 2005). Sev- 
eral attempts to determine it with SINFONI integral-field 
Ha spectroscopy have been performed. For example, Cresci 
et al. (2009) found at z ~ 2.2 a slope consitent with the Bell 



& de Jong (2001) z ~ relationship, while in Vergani et al. 
(2012) the scatter is higher at z ~ 1.5, especially for those 
galaxies with lower rotational velocities. 

We attempt to shed some light into the TFR for our 
sample. The number of objects is not high, and their mor- 
phological nature is not perfectly constrained, but in spite of 
these limitations this is one of the few cases where this rela- 
tionship could be obtained for IFU observations of massive 
galaxies at high redshift. Moreover, it helps understanding 
how the mass and the rest of kinematic properties of our 
sample are connected. 

We show the K s band TFR for our sample of galaxies 
in Fig. [7] using the maximum rotational velocity retrieved 
in our modelling. As we did in the previous plots, we add 
SINFONI Ha measurements for the massive galaxies in the 
SINS survey (Forster-Schreiber et al. 2009, with a detailed 
TFR study in the aforemetioned Cresci et al. 2009). The 
solid line accounts for the local (z ~ 0.2 — 0.3) K s band 
POWIR/DEEP2 relationship found in Fernandez-Lorenzo 
et al. (2010), with the crosses the objects studied in that ar- 
ticle. Overall, our results are similar to Epinat et al. (2009) 
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Figure 7. Tully-Fisher relation with K s absolute luminosity for our massive galaxies at high redshift. The number of each galaxy is 
plotted, to better understand the properties of each galaxy throughout the paper. Symbols have the same meaning as in previous plots, 
with green squares representing the massive galaxies in the SINS sample (Forster-Schreiber et al. 2009). The solid line comes from 
Fernandez-Lorenzo et al. (2010), and it is the local (z ~ 0.2 — 0.3) Tully-Fisher relation derived for POWIR K a band galaxies (which the 
small crosses). This relationship was inferred by weighting the importance of every point by its errors. Galaxies displaying the lowest 
rotational velocity values are the ones which depart more from the local relation. 



and Vergani et al. (2012), as the galaxies which display the 
lowest rotational velocities (POWIR2 and POWIR11) are 
the ones further from the fiducial relation. When plotting 
the results for POWIR4 and POWIR10, they are placed in a 
similar locus, but the reason for this would be their allegedly 
non-massive nature . Whereas for POWIR2 this constitutes 
a further indication of its departure from a pure disk system, 
the interpretation is not so obvious for POWIR11, when 
looking also at the anisotropy plot (Fig. We attribute 
this to the fact that this is the object with the lowest incli- 
nation in our sample, and subsequently it is more difficult 
to constrain this parameter which affects the rotational ve- 
locity determination. 

The stellar mass or baryonic TFR (using stellar mass 
instead of luminosity) has been claimed to be a better proxy 
for studying galaxy mass assembly than the usual TFR. We 
show the baryonic TFR for our sample in Fig. [8] The solid 
line is the local relationship found by Bell & de Jong (2001), 
corrected to our Chabrier IMF, while the dashed line is the 
derived stellar mass TFR at z ~ 2.2 in Cresci et al. (2009). 
We also add the disk galaxies from Conselice et al. (2005), 



separating their sample between z ^ 0.7 and z > 0.7 as 
they did. There is a close resemblance with Fig. [7] We can 
see that all the massive galaxies in the present study occupy 
similar loci, giving further evidence for our disk assumption 
in the galaxy modelling. 



In order to disentangle better the disordered motions of 
the gas we follow the prescriptions in Kassin et al. (2007) 
- see also Cresci et al. (2009) and Vergani et al. (2012) -, 
where they used the parameter So. 5 = yj (0.5 * Vmax) + & 2 , 
arguing that the scatter in the TFR is tighter when tak- 
ing into account the contribution in the velocity dispersion. 
We show this relation in Fig. [9j plotting the highest red- 
shift (0.925 < z < 1.2) relation inferred in Kassin et al. 
(2007). There does not appear to be any conspicuous dif- 
ference from the rest of TFR plots, suggesting again that 
the massive galaxies with the lowest values of maximum ro- 
tational velocity have properties that are more difficult to 
match with the assumption of a pure disk-like nature. 
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4 DISCUSSION & CONCLUSIONS 

We present a SINFONI study of a sample of massive galaxies 
(Mgteiiar ^ 10 11 h^g Mq) at z ~ 1.4 selected by their stellar 
mass and emission lines (with EW[OII] > 15 A), in order 
to understand the kinematics/secular motions of this galaxy 
population and to constrain their rotational nature. This is 
a matter of debate after recent photometrical studies (e.g. 
van der Wei et al. 2011, Buitrago et al. 2013 among oth- 
ers) found that high redshift massive galaxies are predomi- 
nantly disk like. There are also Ha detections using slitless 
spectroscopy (cf. 3D-HST Survey; van Dokkum et al. 2011) 
showing intense star formation in massive high-z galaxies. 
Our work is an attempt to clarify the diversity of properties 
within these galaxies, and whether they are better described 
kinematically by a disk-like or a spheroidal-like population. 

In this study we have carefully chosen 10 massive galax- 
ies with available deep Keck spectroscopy and K s band imag- 
ing from the POWIR/DEEP2 survey. VLT/SINFONI H- 
band observations with a very good (average 0.56 arcsec) 
seeing, enabled us to build Ha kinematic maps. We fit ro- 
tating disk models to their velocity fields that allowed us to 
derive rotation velocities and correct the velocity dispersion 
maps from beam smearing. Hence we minimize potential er- 
ror sources as the broadening of the spectral lines by velocity 
shear. 

An open question is whether we could generalize our 
conclusions to the whole massive galaxy population at 
z~1.4. To answer this we need first to ascertain whether our 
sample is unusual for showing [EW]> 15 A, and thus non 
negligible star formation rates, as both Ha and [Oil] come 
from ionized gas by star formation in HII regions. Among 
studies similar to ours, Twite et al. (2012) analyse Ha star 
formation in M steUar >10 10,5 M Q galaxies at 0.4 < z < 1.4. 
Although their sample was not complete in stellar mass, they 
found a significant drop at lower redshift in the number of 
massive galaxies with detected Ha emission. 

Although they do not sample the same redshift range 
as our study, Sobral et al. (2011) find that the the frac- 
tion of massive galaxies with M s teiiar ^ 10 11 Mq detected 
in Ha is ~ 15% at z — 0.84. These systems have equiva- 
lent widths greater than ~15 A, which typically translates 
into star formation rates of 5-10 Mqyt~ 1 (P. Best private 
communication). Finally, Fumagalli et al. (2012), working 
with slitless spectroscopy in the 3D-HST survey, have found 
(based on their Figure 2) that at z~1.4 massive galaxies 
typically show equivalent widths similar to the ones in our 
sample. Summarizing, despite (on average) a level of star 
formation and thus a certain amount of Ha emission is ex- 
pected, our conclusions should be taken carefully as it is 
not perfectly clear yet whether all massive galaxies have the 
same characteristics of our sample at this specific redshift 
(z-1.4). 

The massive galaxies in our sample show remarkably 
ordered rotational velocity gradients and high velocity dis- 
persions, as has been seen before in star-forming but less- 
massive 3D spectroscopy samples at high redshift. All galax- 
ies from our sample show Vm. ax /a > 1, where this ratio is 
on average 3—4 in most cases (Figure [4| . The velocity dis- 
persions, although smaller than the typical values reported 
for the stellar component in massive galaxies at these red- 
shifts (see for example Cenarro & Trujillo 2009), resemble 



those found for ionised gas kinematics. These large V m ax/<J 
ratios are however at odds with local Universe counterparts, 
which either display Vmax/c < 1 (e.g. Emsellem et al. 2011) 
for early-type galaxies or V m ax/<J ~ 6 — 14 (e.g. Epinat et 
al. 2010) in case of IFU studies of Ha emission in spirals. 

We agree with the explanation given in previous high 
redshift 3D spectroscopy studies (Law et al. 2009, Forster- 
Schreiber et al. 2009, Epinat et al. 2012) such that, at high 
redshift, galaxy formation and evolution is a more turbu- 
lent process because of the larger amounts of cold gas in- 
volved, which at the same time leads to higher star forma- 
tion rates than in the present day Universe. Likewise, we 
observe that major merging is indeed occurring in our sam- 
ple (see POWIR4 or POWIR10 galaxies). However most of 
the gas should be accreted either via minor merging, whose 
hints are found in multiple galaxies of our sample, or cold 
gas flows along cosmic web filaments (e.g. Conselice et al. 
2013). 

The main difference between our sample and previ- 
ous published datasets (such us MASSIV, LSD/ AMAZE 
or SINS surveys) lays in our high stellar mass selection 
(M a teiiar 10 11 K^q Mq ) . Observationally, we find that our 
sample consists of quite regular velocity fields showing high 
rotation. As stated in Epinat et al. (2012) when discussing 
their disk galaxies, this fact implies that the most massive 
disks seem to be stable objects even at early cosmic times. 
We have compared our massive galaxy sample with other 
galaxies in the literature selected by mass. Strikingly, the 
conclusions remain the same. In addition, we have noticed 
that less-massive galaxies (Law et al. 2007, Wright et al. 
2007, Epinat et al. 2009, Forster-Schreiber et al. 2009, 2011, 
Epinat et al. 2012) contain a higher percentage of clumpy 
or distorted objects than our sample. 

We conclude that massive galaxies acquire more rapidly 
a given morphology and gravitational equilibrium than less- 
massive objects, accounting for a type of morphological 
downsizing. We attribute this characteristic to the fact that, 
as we are dealing with the most massive galaxies at these 
redshifts, their high masses help protect them from being 
perturbed. Stellar mass has not only a profound impact in 
the galactic assembly and the star formation history of these 
galaxies, but it is also crucial for understanding their even- 
tual morphological development, whereby they progressively 
obtain the observational properties of the massive galaxies 
we find in the nearby Universe. 

Future NIR high resolution photometry over larger sam- 
ples of massive galaxies shall contribute to corroborate this 
scenario. These objects should be the base of surveys taking 
advantage of new generation multi-IFU spectrographs which 
will increase by a high factor the number of galaxies with 
available kinematic information. 
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Table 2. Modelled kinematic parameters of our sample 



Name 


Velocity dispersion 


Max. rotational velocity 


V max /a 


Vel. disp. in r e 


Max. rot. vel. in r e 


(V m „/a) e 


Inclination 


Classification 




km s -1 


km s -1 




km s 


km s — 1 




° (degrees) 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


POWIRl 


70 ±41 


236 ± 19 


3.37 ±0.29 


71 ±40 


236 ± 15 


3.32 ±0.21 


80 


D 


POWIR2 


72 ±30 


117 ± 16 


1.63 ±0.24 


89 ±7 


117 ±9 


1.31 ±0.09 


64 


D 


POWIR3 


42 ±26 


170 ± 25 


3.96 ±0.62 


18 ±29 


76 ± 19 


4.27 ±1.14 


65 


I 


POWIR4 


71 ±30 


95 ±35 


1.34 ±0.50 


74 ±22 


95 ± 13 


1.28 ±0.17 




I 


POWIR5 


131 ± 60 


313 ± 28 


2.40 ±0.22 


172 ± 66 


313 ±37 


1.82 ±0.21 


48 


P 


POWIR6 


43 ±35 


214 ± 29 


4.88 ±0.69 


38 ±32 


214 ± 28 


5.63 ±0.75 


57 


D 


POWIR7 


20 ±23 


141 ± 27 


6.88 ± 1.42 


17 ±22 


141 ± 26 


8.16 ±1.58 


63 


P 


POWIR8 


60 ±23 


278 ± 27 


4.60 ± 0.47 


87 ± 18 


99 ±21 


1.14 ±0.24 


61 


D 


POWIR10 


59 ±26 


113 ± 18 


1.92 ±0.32 


53 ±22 


114 ± 10 


2.13 ±0.19 




I 


POWIR11 


45 ±23 


111 ± 15 


2.43 ±0.36 


24 ±22 


111 ± 10 


4.53 ±0.46 


31 


D 



Notes. (1) Name of the galaxy (2) ' 1/error 2 '-weighted average velocity dispersion from the modelled velocity dispersion after removing 
the beam smearing (3) Maximum rotational velocity from our rotational velocity modelling (4) Ratio between the maximum rotational 
velocity and velocity dispersion (5) ' 1/error 2 '-weighted average velocity dispersion from the modelled velocity dispersion within one 
effective radius after removing the beam smearing (6) Maximum rotational velocity at one effective radius from our rotational velocity 
modelling (7) Ratio between the maximum rotational velocity and velocity dispersion at one effective radius (8) Inclination as measured 
from GALFIT analysis. Note that for POWIR4 and POWIR10 values are fixed to 60°, as for these two cases the observed objects are 
thought not to be the targeted massive galaxies (9) Final kinematic classification for our massive galaxies: D for rotating Disks, I for 
Interacting galaxies and P for Perturbed rotators. 

Table 3. Radii used in our calculations 



Name 

(i) 


Model radius 
(2) 


Model radius 
kpc 

(3) 


Ha maximum radius 
(4) 


Ha maximum radius 
kpc 

(5) 


Effective radius 
kpc 
(6) 


Pixels taken 
(7) 


Radius taken 
kpc 
(8) 


POWIRl 


0.12 


1.05 


0.62 


5.26 


3.36 


3.5 


3.68 


POWIR2 


0.12 


1.05 


1.10 


9.27 


2.59 


2.5 


2.63 


POWIR3 


0.70 


5.90 


1.22 


10.30 


3.06 


2.5 


2.63 


POWIR4 


0.12 


1.05 


1.05 


8.83 


4.18 


3.5 


3.69 


POWIR5 


0.43 


3.64 


1.03 


8.68 


3.77 


3.5 


3.68 


POWIR6 


0.16 


1.35 


1.44 


12.10 


4.81 


4.5 


4.74 


POWIR7 


0.12 


1.05 


0.88 


7.41 


3.85 


3.5 


3.69 


POWIR8 


0.88 


7.43 


1.44 


12.17 


2.58 


2.5 


2.63 


POWIR10 


0.12 


1.05 


1.06 


8.90 


2.87 


2.5 


2.63 


POWIR11 


0.12 


1.05 


0.97 


8.19 


2.83 


2.5 


2.63 



Notes. (1) Name of the galaxy (2) Radius of our kinematic model in kpc (3) Radius of our kinematic model in arcsec (4) H« maximum 
extent in arcsec (5) Ha maximum extent in kpc (6) Effective radius as calculated by the disk-like relation for massive galaxies in Table 
2 of Buitrago et al. (2008) (7) Pixels taken as effective radius according to previous column (note that the 0.5 is added as we start from 
the kinematic center). (8) Equivalent in kpc of the previous column. 



APPENDIX A: INDIVIDUAL GALAXY 
OBSERVATIONS 

In this appendix we present a detailed description of each 
massive galaxy within our survey and their observed kine- 
matics. Before going into this, we preface this discussion by 
showing the model maps in Fi eures [A 1 1 and [A2l Each galaxy 
explanation compares the various maps (cf. Figures IA3I - 
|AT3|) , namely the parent POWIR/DEEP2 K 3 band imag- 
ing, the Ha flux images from SINFONI, H-band continuum, 
SNR map, radial velocity (and its residual), observed ve- 
locity dispersion (and the inferred velocity dispersion after 
removing the beam smearing). 

In our maps, north is up and east is left. Ha contours 
are overplotted in all of the maps (with decrements of 10% 
in flux between adjacent contours), in order to facilitate the 
reader in determining which spaxels belong to the galaxies. 
The axes show sizes both in kpc and arcsec. The kinematic 
centers used for our models are located in the spaxels with 
the maximum flux in the continuum maps, and are high- 
lighted by a cross. Seeing values are presented by circles with 
their FWHM corresponding size in the left bottom corner 
of the SINFONI maps. In the kinematic maps, the coloured 



spaxels shown are all above a certain threshold in SNR which 
is written in Table [T] Column (S/N) maps . However, this 
threshold was not applied for the continuum maps, in or- 
der to understand the origin of the Ha emission. Finally, we 
attempted to quantify the existence of AGN sources within 
our sample by means of obtaining the N2 indicator - the 
logarithm of the [Nil] /Ha ratio; see Queyrel et al. (2009) 
for details -. For POWIRl, 3 and 7 this could not be ac- 
complished as the A6583A [Nil] line is located over OH sky 
lines (see Table Q. We proceed to describe any interesting 
feature of our sample in the following galaxy subsections. 

Al POWIRl 

The maps belonging to this galaxy are shown in Figure IA3I 
This galaxy, albeit a clear detection, is a very compact sys- 
tem both in the K s band and in the Ha image. The explana- 
tion for observing only few spaxels is its inclination, which 
is the highest of the sample. When looking at the contin- 
uum image other features appear. There is a flux stripe in 
the right side which is spurious as it extends from the left to 
the right part of the detector, having no counterpart in both 
K s and Ha maps. However, we notice two blobs which seem 
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Table 4. Masses inferred for our sample & N2 calibrator 



Name Stellar mass Dynamical mass Mass enclosed in Ha max. radius Asymmetric drift correction N2 calibrator Notes about N2 calibrator 

10 10 hf 2 M e 10 10 h7„ 2 M s 10 10 h- 2 M o 10 10 fef 2 M o 

(1) (2) (3) (4) (5) (6) (7) 



P0WIR1 


17.1 


9.7 


6.9 


2.8 




OH line over [JV//]6584Awavelength 


P0WIR2 


10.4 


30.6 


3.0 


27.6 


-0.31 




P0WIR3 


14.3 


16.6 


6.9 


9.7 




OH line over [jV-Z7]6584Awavelength 


P0W1R4 


25.7 


10.9 


1.9 


9.0 


-0.59 




P0WIR5 


21.2 


54.9 


19.8 


35.0 


-0.22 




P0W1R6 


33.6 


19.6 


13.0 


6.6 


-0.31 


OH sky line residual increases this ratio 


P0W1R7 


22.1 


3.9 


3.4 


0.5 




OH line over [JV//]6584Awavelength 


P0W1R8 


10.4 


66.3 


22.(1 


44.3 


-0.60 




POWIR10 


12.7 


16.0 


2.7 


13.3 


-0.33 


The value for the massive galaxy is N2 = —0.15 


P0WIR11 


12.4 


8.8 


2.4 


6.4 


-0.38 





Notes. (1) Name of the galaxy (2) Stellar mass from the parent POWIR/DEEP2 survey (3) Dynamical mass, as calculated in Section 
13.31 (4) Enclosed mass term, as calculated in Section 13.31 (5) Asymmetric drift correction, as calculated in Section 13.31 (6) N2 calibrator 
as in Queyrel et al. (2009) (7) Observational notes 



to be real. They have an angular size comparable with the 
seeing of this observation, which reassures us of their detec- 
tions. The fainter one, near the northern edge of the primary 
galaxy, may help us understand why this galaxy shows large 
velocity dispersion values close to the blob, as this may sig- 
nal a minor merging event. Conversely, the brightest spot on 
the south west of the continuum image of the galaxy is not 
associated with any Ha emission, and it does not cause any 
significant distortion of the main galaxy. Kinematic models 
show a regular rotational gradient and a fairly high (~ 70 
km/s) 1/error 2 velocity dispersion. Although it is tentative 
to identify this galaxy as a merger, we prefer to classify it 
as very inclined disk galaxy, because of its large rotational 
velocity and ordered velocity field without any substantial 
disruption. 

A2 POWIR2 

The maps belonging to this galaxy are shown in Figure lA4l 
The Ha flux map covers in this case the whole K s image of 
the galaxy. The continuum center and Ha center are well 
aligned. The Ha map peaks in the center, and there is also 
a very bright group of spaxels on the north west part of the 
galaxy. This sharp feature is most probably caused by a cos- 
mic ray as it has no counterpart both in the continuum and 
SNR images. The galaxy looks asymmetrical in its southern 
part, where we found an enhancement of the SNR as well. 
We are perhaps witnessing an accretion of a minor object. 
But what we can state is that this galaxy looks like a relaxed 
and ordered system. The Ha line lies in a spectral region far 
away from any sky line and thus the [Nil] line is clearly 
identified, giving a ratio of [NII]/Ha = 0.49. This 1:2 ratio 
between Ha and [Nil] is preserved within the external spax- 
els, and may indicate the presence of an AGN. We consider 
this system as an early disk-like galaxy or spheroid with a 
disk, as its large velocity dispersion and low V ma x/<J ratio 
reveal. Nevertheless, the interpretation of a disk galaxy is 
favoured at the light of the strong velocity gradient. 

A3 POWIR3 

The maps belonging to this galaxy are shown in Figure lA"5l 
Both K s and continuum maps show an elongated structure 
with a diagonal shape from east to west. Hacountours do 
not exactly overlap with the galaxy continuum, and none of 
the brightest Ha spots coincide with it. The Ha center and 



the continuum center are located in different places. Regard- 
ing the continuum, we rely on its location, but not on its 
exact shape as it may be affected by the aforementioned de- 
tector problems. This is a clear case of a disturbed object, 
but it is remarkable that even in this case, the rotational 
field is quite clear. The lack of any neighbouring galaxy and 
stretched shape are evidence for an ongoing merger. With 
the data we have we cannot add anything to this discus- 
sion. It is important to stress that this is one of the galaxies 
which was observed for half of the integration time, and for 
this reason it does not have a well-defined shape. Taking 
all the available probes into account we define it as dis- 
turbed/merging galaxy. 

A4 POWIR4 

The maps belonging to this galaxy are shown in Figure lA"6l 
If one looks either at the K s band map or at the contin- 
uum map a bright galaxy appears, and by its eastern side 
a very elongated arc-shaped feature of Ha , which has some 
very weak continuum as well. In this case, it seems that 
the emission comes from a minor object whose gas has been 
stripped or conversely a fan of stars coming from the main 
object. The Palomar image also shows this feature, indicat- 
ing that the merging interpretation is favoured. Our kine- 
matics are thus inferred for the Ha visible object which was 
not the primary target of our observations. We find a rota- 
tional velocity field, but it is not large (perhaps due to its 
non-massive nature) and a comparatively big velocity dis- 
persion consistent with its interacting nature. Its morphol- 
ogy is clumpy, with at least two clumps being visible. We 
catalog it as disturbed/merging galaxy. 

A5 POWIR5 

The maps belonging to this galaxy are shown in Figure IA7I 
By looking at the flux images, this galaxy appears as a blob 
with a tail in its upper part. There might be companions 
according to the POWIR/DEEP2 imaging, but we cannot 
state this with certainty. Its Ha emission extends over the 
K s image and the continuum image. The maximum rota- 
tional velocity is very high, as is the velocity dispersion, and 
the result is that the dynamical mass is the second largest 
in our sample. The velocity dispersion must be artificially 
increased by a contribution from a sky line in the redder 
part of the spectrum. The pipeline is meant to remove OH 
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sky lines, but there are sometimes (as in this case) residuals. 
When the Ha line and the sky line are close to each other 
it becomes very hard for our IDL programs to disentangle 
them, increasing somewhat the final result, although not sig- 
nificantly, as our routines were able to resolve the gradient in 
the radial map. We find [NII]/Ha= 0.60, which is puzzling 
as it is a large ratio and we note that the [Nil] is detectable 
all over the galaxy spaxels, and not just concentrated in the 
center, as we would expect for AGN emission. The prop- 
erties of this massive galaxy cannot be perfectly explained 
by a disk-like object and for this reason we classify it as 
perturbed galaxy. 

A6 POWTR6 

The maps belonging to this galaxy are shown in Figure IA8I 
This is the brightest galaxy in our POWIR imaging. The 
Ha appearance of the galaxy is different from either the con- 
tinuum or the K s band image. Remarkably, we do not detect 
Ha emission in the central area, appearing as a hole in the 
Ha flux map and four knots or clumps surrounding it. This 
hole overlaps with the center of the K s and continuum im- 
ages. Similar cases are seen in Epinat et al. (2010), for ex- 
ample after redshifting the local galaxy UGC04820 where its 
ring nature produces the effects observed in our flux maps. 
This is typically found in early- type spiral galaxies. The 
clumps are also conspicuous in the velocity dispersion map. 
Both theoretical expectations and other recent work agree 
with a likely clumpy phase in galaxy formation at high-z 
(Forster-Schreiber et al. 2011b and references therein). We 
identify four of them in the Ha map, which match the ve- 
locity dispersion enhancements. Another piece of evidence 
for a star-forming disk is that it is strongly rotationally sup- 
ported (Vmax/a = 4.88 and V ma x = 214 km s _1 ). Hence, we 
classify it as a disk galaxy. 

A7 POWIR7 

The maps belonging to this galaxy are shown in Figure IA9I 
In this case, one of the images of the galaxy fell in the 
borders of the SINFONI detector. This is the reason for 
its low SNR in all its spaxels. However, K a POWIR/DEEP 
imaging, Ha and continuum maps overlap well, and show 
a slightly distorted and clumpy galaxy. Rotation seems to 
play an important role in the support of this galaxy (with 
Vmax/& = 6.88), but the low SNR prevents us from drawing 
any further conclusions. Because of its irregular shape, we 
consider this is a perturbed galaxy. 

A8 POWIR8 

The maps belonging to this galaxy are shown in Figurc lAlOl 
Our best seeing (0.42 arcsec) observations are for this sys- 
tem where we retrieve a very clear disk. A bulge component 
in the center may be present, as the central region is very 
bright in both Ha and K s imaging, and displays a large ve- 
locity dispersion of a ~ 130 km/s. This agrees with the 
renditions of the anisotropy plot (which will be explained in 
the Section[32}, as this system is found near a locus close to 
the local fast rotators. However, an alternative explanation 
could be that our models were not able to resolve completely 



the rotation in the inner part of the galaxy, in which case 
we would have a ring distribution of gas as in POWIR6. 
The rotational support is very strong {Vmax/o- = 4.60). Its 
dynamical mass is also very large (6.63 x 1O 11 M0) and re- 
liably measured due to its high rotational velocity gradient, 
the identification of which was helped by the excellent see- 
ing of this object's observation. We classify this galaxy as a 
rotating disk. 

A9 POWIR10 

The maps belonging to this galaxy are shown in Figures 
IAlll and lA12l We present two different set of figures for this 
galaxy. Our K s and continuum images show two galaxies. 
Most of the flux in these two maps comes from the object 
in the southern part of the image (Figure lAllf) . although 
Ha emission is mainly found in the galaxy on the north (Fig- 
ure [Xl2}. This system has some similarities with POWIR4, 
specifically in that most of the Ha emission does not come 
from the position where we had expected to find the mas- 
sive galaxy. For this reason, the Ha kinematic parameters 
listed on the tables as POWIR10 are associated with the 
object identified as the non-massive galaxy, as kinematics 
could only be reliably retrieve from it. As a matter of fact, 
we found strong emission at the [Nil] wavelength for the cen- 
tral spaxels of the southern object. This makes us suspect 
that, as it is nearly a point source, this system harbours an 
AGN on its center, and the Ha emission at its sides might 
be related with shocked gas by outflows of material com- 
ing out from it. Regarding the non-massive galaxy within 
our imaging, we are able to detect a hint of continuum and 
very clear kinematics. Due to the low values in the kine- 
matic velocities for this second galaxy we conclude that the 
photometric mass is only derived for the AGN. For the total 
system, due to its distortion/clumpiness, the velocity dis- 
persions enhancement in between the two objects and of 
course since Ha is detected at the same wavelength for the 
two subsystems, we classify this object as a merger. 

A10 POWIR11 

The maps belonging to this galaxy are shown in Figure rAl3l 
This is a clearly detected galaxy which was only observed 
half of the time as the other systems. However, Ha emission 
show an extended galaxy, in agreement with the K s and con- 
tinuum images. It is most probably close to being face-on, as 
its inclination is the smallest (30°) in our sample. The SNR 
decreases in its western side because of the presence of a sky 
line at the same side of the spectrum where the Ha line is 
found. Nevertheless, the Ha line is clearly distinguished in 
all the spaxels up to the faint outskirts, making the rota- 
tional velocity gradient easy to detect. This was the reason 
behind setting S/N = 2 as a threshold for the galaxy maps. 
Overall, this galaxy presents a disturbed discoidal shape and 
a couple of bright Ha spots or clumps in its central part. The 
absence of any merger signature makes us conclude that this 
is a disk, although rather perturbed or turbulent, as indi- 
cated by its fairly low V m ax/<J ratio. 
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Figure Al. Montage with the kinematical models. First row are the rotational velocity maps, and the second the velocity dispersion 
maps. From bottom to top, the columns are POWIR1, POWIR2, POWIR3, POWIR4, and POWIR5. Please note that the full descriptions 
of our maps are written in the Appendix A and in the caption of the first galaxy montage (Figure I A3ft . 
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Figure A2. Montage with the kinematical models for the rest of the massive galaxies not plotted in the previous Figure fAll . First 
row are the rotational velocity maps, and the second the velocity dispersion. From bottom to top, the columns are POWIR6, POWIR7, 
POWIR8, POWIR10 (non-massive galaxy), POWIR10 (massive galaxy) and POWIR11. 
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Figure A3. From left to right, from top to bottom: K s POWIR survey image of the galaxy [K s ], H-band SINFONI continuum image 
[Cont], SINFONI Ha flux image [Ha], signal-to-noise map for every spaxel [SNR], observed radial velocity map [Rot j, s ], residual rota- 
tional velocity map [Rot re3 ] (after subtracting quadratically the rotational velocity model from Figure [ATI in this case, otherwise Figure 
I A2I > . observed velocity dispersion map [ct j, s ] and residual velocity dispersion map [oy es ] (after subtracting quadratically the velocity 
dispersion model from Figure |A"T1 in this case, otherwise Figure lA"2ll , In our maps, north is up and east is left. Ha contours are overlapped 
in all the maps (with decrements of 10% in flux between adjacent contours) in order to facilitate the reader to know which spaxels belong 
to the galaxies. The axes show sizes both in kpc and arcsec. A cross appears in the place we have set our galaxy center (see beginning 
of Section [X] for an explanation) and the seeing is represented in the lower left part of the SINFONI maps by a circle whose size is its 
FHWM. 

POWIRI - Disk-like galaxy. Comments: Clear Ha line in all the spaxels above the signal-to-noise threshold, velocity dispersion enhance- 
ment due to minor merging. 
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Figure A4. POWIR2 - Disk-like galaxy. Comments: Very extended Ho emission, which is very clearly detected. There is an elongation 
and a signal-to-noise enhancement in the southern part, arguably because of a merging episode. The group of spaxels in the top right 
display high velocity dispersion are an artefact, but not the high values in its center. This latter feature points towards the development 
of a spheroidal component. 
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Figure A5. POWIR3 - Merging/Interacting galaxy. Comments: K s and continuum images are vaguely related but not with the H« map. 
Most probably this is an ongoing merger, and this would explain its elongated shape. Rotation is found but it shows complex kinematics, 
as expected given its interacting nature. 
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Figure A6. POWIR4 - Merging/interacting galaxy. Comments: Both K s and continuum maps show a bright object located at the right 
side of the Ha detection. We interpret this as the massive galaxy (which is devoid of Ha emission) is interacting with another object, 
whose gas is being stripped or very perturbed. For the Ha detection, which is the non-massive galaxy, non-ncglectable values of rotation 
and velocity dispersion are retrieved. 
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Figure A7. POWIR5 - Perturbed rotator. Comments: K s , continuum and Ha overlap in the same location, showing a galaxy slightly 
elongated in its top part. This galaxy shows high radial velocity values. Velocity dispersion might be affected by a close sky line, increasing 
slightly its actual value. 
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Figure A8. POWIR6 - Disk-like galaxy. Comments: This ring pattern in the Ha image was observed before in low redshift early-type 
spiral galaxies (see Epinat et al. 2010). We identify four different clumps in Ho and the velocity dispersion map. 



26 Fernando Buitrago et al. 




Figure A9. POWIR7 - Perturbed rotator. Comments: K s , Ha and continuum maps overlap well. This galaxy was not observed half 
of the integration time but was set in the border of the detector, which explains its low signal-to-noise detection. Slightly disturbed 
morphology, with a possible clumpy structure. 



SINFONI/VLT 3D spectroscopy of massive galaxies: Strong rotational support at z ~ 1.4 27 




■15.81 I I I i 

-15.81 -7.91 0.00 7.91 15.81 

X (kpc) 

X (arcsec) 

-1.BS -0.94 0.00 0.94 1.88 



I 

^1 




■15.81 I v ~ t I I I 

-15.81 -7.91 0.00 7.91 15.81 

X (kpc) 



I 





-15.81 -7.91 0.00 7.91 15.81 

X (kpc) 




obs 



■ 



-15.81 -7.91 0.00 7.91 15.81 

X (kpc) 



Rot 



res 



-15.81 -7.91 0.00 7.91 15.81 

X (kpc) 



a 



obs 



-15.81 -7.91 0.00 7.91 15.81 

X (kpc) 



a 



res 




15.81 -7.91 0.00 7.91 15.81 

X (kpc) 



Figure A10. POWIR8 - Disk-like galaxy. Comments: Clear and extended disk in all the images. High radial velocity values, with a 
large velocity dispersion in its center, which is a hint of a bulge component. 
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Figure All. POWIR10 - Merging/interacting object. Comments: Similar case as POWIR4. Here we can see the main object, which 
has another blob in its northern part. As most of the continuum and K s signal come from this main object, we identify it as the massive 
galaxy. It show high [Nil] values in its centers, indicative of its AGN nature. The Ho detections at its sides might be related with gas 
outflows. The kinematics are messy, we cannot infer anything conclusive. 
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Figure A12. POWIR10 - Merging/interacting object. Comments: This is the galaxy which was at the north of the previous Figure 
IA11I object. It is more extended, at least in the Ha map. Its radial and velocity dispersion values are not very high, suggesting it 
is not as massive as its partner galaxy. There is a velocity dispersion enhancement in the southern part, where the two galaxies are 
connected / interacting. 
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Figure A13. POWIR11 - Disk-like object. Comments: Close agreement between K s , continuum and Ha light distributions. This galaxy 
could only be observed for half of the nominal integration time, retrieving a low signal-to-noise per spaxel. This is the galaxy with the 
lowest inclination in our sample. Clear velocity gradient as Ha is clearly detected in all the coloured spaxels. 



